Constraining weak annihilation using semileptonic D decays 
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The recently measured semileptonic Ds decay rate can be used to constrain weak annihilation 
(WA) effects in semileptonic D and B decays. We revisit the theoretical predictions for inclusive 
semileptonic -D(a) decays using a variety of quark mass schemes. The most reliable results are 
obtained if the fits to B decay distributions are used to eliminate the charm quark mass dependence, 
without using any specific charm mass scheme. Our fit to the available data shows that WA is smaller 
than commonly assumed. There is no indication that the WA octet contribution (which is better 
constrained than the singlet contribution) dominates. The results constrain an important source of 
uncertainty in the extraction of \Vub\ from inclusive semileptonic B decays. 



An intriguing hint of a possible conflict in the B- 
factory data, which may be a sign of physics beyond the 
standard model, is the roughly 2a difference between the 
value of I Vub \ obtained from a global fit to the CKM pa- 
rameters using sin 2/3 extracted from the CP asymmetry 
in _B — !• Ks and related modes, and \Vub\ measured di- 
rectly from inclusive semileptonic B — Xutv decays, cal- 
culable using an operator product expansion (OPE) 
^. At order Aq^p/rnj^ in the OPE, four-quark operators, 
the so-called weak annihilation (WA) terms, give a signif- 
icant contribution in a phase space region which affects 
all inclusive \Vub\ measurements to some extent. Hence, 
a reliable estimate of the WA contribution is necessary 
to determine whether there is, in fact, a conflict with 
standard model predictions. 

The WA contribution to the total B X^tv decay 
rate |^, ^, and to the charged lepton (or neutrino) en- 
ergy spectrum are calculable in terms of the matrix 
elements of local four-quark operators. However, there is 
so far no first-principles derivation of the WA contribu- 
tion to the double or triple differential spectra. For this 
reason, the extraction of the WA contribution from the 
differential B Xuiv spectra is model dependent, and 
any model-independent bound on the magnitude of the 
WA matrix elements is important. It was pointed out 
by Voloshin |^ that the same matrix elements that enter 
B — >■ Xutv decay can be constrained by the semileptonic 
rate difference between D and Ds mesons, since the B 
and D matrix elements are related by heavy quark sym- 
metry. In this paper we revist the theoretical calculations 
of semileptonic D decays and extract bounds on the WA 
contribution to D and B decays. 

At order Aqq^/tti^, there are dimension-6 four-quark 
operators in the OPE for the semileptonic D decay rate, 
the WA operators. 
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where q = s,d and Cy is the heavy quark effective theory 
charm quark field. The matrix elements B^2^ of these 



operators are defined by 
(Aloi^l^lA) 
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where i = u,d,s labels the fiavor of the light quark in 
the D meson. Compared to the dimension-3, 5, and 
other dimension-6 operators, the matrix elements of these 
operators are enhanced by IGtt^, and contribute to the 
semileptonic decay widths of the three D mesons as 
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relative to the semileptonic width at lowest order in the 
OPE and at tree- level in as , 
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If one assumes factorization and the vacuum saturation 
approximation, then the WA contribution vanishes, since 
Bi — B2 — 1 or depending on whether or not q is the 
same as i. Deviations from the factorization ansatz are 
usually estimated at the 10% level 

The above analysis also holds for B decays, with the 
replacement of D meson quantities by the corresponding 
B meson ones. The matrix elements in the c and b sectors 
are related by heavy quark symmetry. The same matrix 
element estimate for B decays (i.e., \Bi — B2\ = 0.1), 
along with /b ~ 200 MeV Q , implies that the four-quark 
operators contribute ~ 3% to the total B — > XuiD rate, 
making it difficult to accurately determine the WA con- 
tribution from B decays. However, the WA contribution 
to D decay is formally (rrib/rncY enhanced relative to B 
decay, and is comparable to the leading order decay rate. 
To, due to the IGtt^ enhancement. Thus studying WA 
effects in D decay is a good way to constrain the matrix 
elements of the four-quark operators. Note that to de- 
termine the WA contribution to B decays at the 1% level 
only requires the four-quark matrix elements to ~ 30% 
accuracy. Even if this 1 / mj? contribution to D decays is 
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comparable to the leading order rate, this does not neces- 
sarily mean that the 1 /rric expansion breaks down, since 
the WA contribution is the only IGtt^ enhanced contri- 
bution at 0{l/ml). 

It has long been known that the difference in the 
semileptonic branching ratios M nffl 



B{D+ 



Xe+v) = (6.49 ± 0.09 ± 0.11) % , 
Xe+v) = (16.13 ± 0.10 ± 0.29) % , 



(5) 



is mainly due to the lifetime difference, and that the 
semileptonic widths V{D'^ Xe+u) k. T{D+ Xe^ 
are equal to within 3 %. Recently CLEO-c measured 
the Dg semileptonic branching ratio 

B{Dt Xe+u) = (6.52 ± 0.39 ± 0.15) % . (6) 

The expressions for semileptonic -D(5) — > Xtu decays 
are well known in the literature. Schematically, 



TsL — Tec + FwA 
where Fwa is defined in Eq. (p|), and 
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where Ai_2 and pi^2 are the matrix elements of dimension- 
5 and 6 two-quark operators, 7i,2,3,4 are the matrix ele- 
ments of time-ordered products, and r = m'^/m'^. These 
may all be determined from fits to various B decay spec- 
tra []TT|-p^. The complete expression including as cor- 
rections is complicated. In our analysis, we include the 
radiative corrections to order |]l4| in the 1/m^ terms, 
and power corrections to order Aq^p/ml? |1^, |l^. For 
the leading term in the OPE, we include the effect of a 
nonzero strange quark mass, since it affects the semilep- 
tonic width by ~ 6%, and set rUg — )• elsewhere. For 
nonzero r, the pi contribution to Fcg has a logr diver- 
gence as r — > 0. In the OPE with r = 0, this divergence 
is effectively absorbed into the matrix elements of WA 
operators. Including the WA contribution converts the 

spectrum into a plus distribution [|l^ , which integrates 
to zero and gives the 77pi/(6m^) contribution to the to- 
tal semileptonic rate in Eq. (^). 

The terms in Fee are all independent of the flavor of the 
spectator quark in the heavy meson, and so give equal 
contributions to Fsl for all three D (or B) mesons. The 
leading term in Fsl which depends on the flavor of the 
spectator quark, and thus produces a difference in the 
semileptonic partial widths, is Fwa- Fwa depends on two 
independent matrix elements in the fiavor SU{3) limit, 
since the operators in Eq. (P have an SU (3) singlet and 
octet part, each of which yield one invariant with the 



two D fields. We can write the decay rates of the three 
D mesons in terms of these two parameters as 



SL — Fcc + Go^|Vc; 

F^L*^ =rce + Go(|K. 
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Here we normalized the weak annihilation contribution 
to the observed semileptonic decay rate. 



Go = 0.16pS-l =Fo(TOrcf) 
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where mref = 1.357 GeV and oq and as are dimensionless 
numbers proportional to the singlet and octet matrix el- 
ements, 



B. 



B 



= as + -^(ao - as) 



Gp mj ff) mp 
12ttGo 

The size of ao,8 is then (approximately) the fraction of 
the meson semileptonic width due to WA. 



The difference F, 
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SL ~ ^ SL suppressed by 

\Vcd/Vcs\'^. Neglecting this correction, it is straightfor- 
ward to extract as from the measured difference of the 
semileptonic widths of the D^'^ and the Dg, as proposed 
in Refs. [|[ |l^. This difference combined with any of 
the individual semileptonic widths also allows ao to be 
extracted, but this requires a reliable computation of 
Fee fisl] . Since the charm mass is not particularly large 
compared with nonperturbative QCD scales, this com- 
putation suffers from both large perturbative and 1/m" 
corrections, limiting the precision with which WA can be 
studied in charm decays. 

The leading perturbative corrections to the semilep- 
tonic decay widths are given by a perturbation series mul- 
tiplying the free-quark decay width Fo given in Eq. (Q). 
This perturbation series depends on the choice of scheme 
for TOe, and rric could be determined from the B decay 
data using the method of Ref. [Q. However, it is well- 
known that the perturbation series for this leading term 
in the OPE is badly behaved when the rate is expressed 
in terms of the charm quark pole or MS masses. In- 
cluding the known results up to order [|4| and using 
as — 0.35 gives the series 



Fo [mr'°] " 
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Fo[mc(TOe)] 



1 - 0.269 e - 0.360 cIlm + 0.069 



= 1 + 0.474 e + 0.513 Eblm - 0.142 - 



(12) 
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respectively. Here e = 1 counts the order in the perturba- 
tion series, and the BLM subscript refers to the e"/3g~^ 
terms in the perturbation series. [In Eqs. (n3) - (n6h we 
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set rUs — i> for simplicity; this has no effect on our dis- 
cussion.] These series are poorly behaved, and do not 
appear to converge. 

The bad behavior of these perturbation series is un- 
derstood theoretically from b decays, and arises from a 
poor choice for the heavy quark mass. A better behaved 
series is obtained by using a threshold mass scheme, such 
as the 15 ||3-||3|, kinetic H] or PS H mass schemes. 
As observed already in |^ , the perturbation series relat- 
ing Tec to the 15 mass is reasonably well-behaved, and 
extracting rric using the method of Ref. pO| 

1 - 0.133 e- 0.006 e|LM - 0-017 e^. (14) 



1.2. Radiative corrections to the l/m^ terms are 



The series is less well-behaved in the PS or kinetic 
schemes (defining both with a 1 GcV factorization scale). 
For the PS scheme we find 

= 1 + 0.262 e + 0.217 m - 0.079 e^, 



ro[mPS(iGeV)] 



ro[mki"(lGeV)] 



(15) 

while for the kinetic scheme, as previously noted p6| , the 
series is considerably worse 

1 + 0.628 e + 0.631 el^^ ~ 0-126 e^. 

(16) 

In addition to the uncertainties in the above series, 
there will be additional uncertainties in extracting a 
charm quark threshold mass from other physical quanti- 
ties, such as moments of B decay spectra ||ll|, |l^. Since 
the charm quark mass is an intermediate quantity which 
is not required for our analysis, we can minimize this 
source of theoretical uncertainty by bypassing any choice 
of charm mass scheme, and instead directly relate the 
semileptonic D decay widths to the values of and 
A = rrib — rric extracted from a global fit to B decay 
spectra. From Eq. (12) and the relation between the h 
quark pole and 15 masses [p2| 



m 15' 



= 1 + 0.011 e + 0.019 e^LM - 0.003 + . . . , (17) 



we find the reasonably well-behaved perturbation series 



= l-0.075e-0.013e|LM-0-021e^ (15 



FoK' 



using 



4.7 GeV, A = 3.4 GeV, a,(mf,) = 0.22, 
and, as in the previous expressions, we have continued 
to set nis to zero. We will therefore use this method to 
determine the D semileptonic widths theoretically. 

The masses 'm\^ and mt, — nic and HQET parameters 
•^1,2, Pi, 2 and 71,2,3,4, as well as their correlated uncer- 
tainties, are obtained using a fit to the B decay spec- 
tra ||ll[ |l^. The values for D decay are related to those 
for B decay by renormalization group evolution between 
rub and rric- Xi is not renormalized due to reparametriza- 
tion invariance p7[|, while A2(toc) = '*cA2(to6), with 



computed in Refs. |2S|-|30|1 . Since they are small and were 
not included in the B decay fits, we neglect them here. 

The errors from the fits include the experimental un- 
certainties, as well as additional theoretical uncertain- 
ties due to neglected higher order corrections, as given in 
Ref. ||l^, We treat the B and D decay calculations 
as independent. Thus the B decay fit results will be held 
fixed (at order Eblm) while we vary the order of the D 
decay results between tree-level and e^. While this may 
be formally inconsistent, numerically, the as and l/mq 
corrections are significantly larger for D than for B decay. 

Using the value of Fee obtained as discussed above, and 
fitting to the experimentally measured rates in Eqs. (|^) 
and (|6|) gives the WA annihilation parameters oq and ag. 
Figure □ shows the 90% CL contours at tree level, order 
e, and order e^. The best fit parameters at order are 

ao = 1.25 ±0.15, 

as = -0.20 ±0.12, (19) 

where the error is from the order fit. The series of a" 
corrections to Tec are flavor independent, and lead to a 
shift in flQ depending on the order in e, but do not af- 



fect ag, which can be determined from F 
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^ cc 

cancels in this difference, so ag is not affected by the con- 
vergence of the as expansion. The shift in ao between e 
and is 0.06, which is smaller than other uncertainties. 
The expansion in 1 /nic is also not as rapidly convergent 
as in the B meson system, so there are significant uncer- 
tainties which mainly affect ag. We find that the 1/m^ 
and 1/m^ terms in Eq. (0) contribute roughly —50 % and 
+35 % to the semileptonic widths. These corrections are 
much larger than the corresponding ones for the hadron 
masses, because of the larger coefficients of A2, pi, and 
P2 ■ One could estimate the uncertainty corresponding to 
these large corrections by including an additional error 
of 0.2 in ao, which is half the 1/m^ term. In the Nc ^ 00 
limit, the meson sector of QCD has a U{3)q ® U{3)q 
symmetry | |3l|] and this implies that ao — ag (see, e.g., 
Ref. [Q), which is shown as the black line in Fig. |^. 

Neglecting Cabibbo-suppressed terms, the correspon- 
dence between our notation and that of Ref. [01 is 



B[" 



(20) 



ref 



and the same equation with ag — > ao and the non- 
singlet Bi 2 replaced by the singlet ones. Taking fn w 
200 MeV [| gives (Bl'" ~ B'^^) « -a8/4.3 « 0.05 ± 0.03, 
which is somewhat smaller than (although consistent 
with) the simple estimate (B]"^ - B^'^) ~ 0.1 in §, |l|l. 

The linear combinations (ao ± 2ag)/3 and (ao — ag)/3 
that contribute to the decay rates in Eq. (|9|) are 

(ao + 2ag)/3== 0.29 ±0.10, 
(ao-ag)/3 = 0.48 ±0.06, (21) 

where only the fit uncertainty is quoted, as discussed 
above. The 90% confidence level contours in these vari- 
ables are shown in Fig. |l|. While there are significant 
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FIG. 1. The 90% CL contours for fits at order e'' (dotted green), (dashed blue), and (solid red). Note that as is not 
affected by the order in e. The thin dot-dashed black line is the large- A^c relation, ao = ag. 



uncertainties in the fit result for the WA contribution in 
Eqs. (^9|) and (^), it still has important implications for 
B and D decays and the determination of \Vub\- 

It has often been assumed that the WA term where the 
light quark in the operator matches that in the heavy me- 
son is much larger than when the light quarks differ, i.e., 
|ao + 2as\ ^ \ao — as\- Indeed, the central values of our 
results suggest that the WA contribution to i?° decay is 
larger than that to decay. The WA matrix element in 
which the light quark field of the operator is contracted 
with the spectator quark in the heavy meson is helic- 
ity suppressed by rnj/rn^, where rm is the lepton mass, 
and gives a contribution of relative order AQQj-,m^/m^ to 
the decay width. Other diagrams, in which the specta- 
tor quark is not annihilated by the four-quark operator, 
are of relative order AQ(-;p/m|?. In a quark model, they 
would contain additional suppression factors from gluon 
exchange to connect the spectator light quark with the 
rest of the diagram, but nothing as small as m^/m^. 

The D meson lifetimes also depend on the WA matrix 
elements through both the semileptonic and non-leptonic 
decay rates. The non-leptonic rates depend on two addi- 
tional color octet operators, and the behavior of the as 
perturbation series is even worse than for the semilep- 
tonic case. Neglecting the color octet matrix elements 
and SU{3) violation (as before), one would predict 
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SL 



total total 



;C+C_ cos2 6'c 



0.3, 



(22) 



where C_ = C^^ = [as{mc) /as{mw)]^^^^^ , and we have 
used C- = 1.6 and C+ — 0.8 for the numerical values. 
The D branching ratios Eqs. (||) and (||) and the lifetimes 
yield 0.07 ± 0.02. This shows that there must be some 
other large contribution to the nonleptonic decay rates, 
e.g., large color octet matrix elements, corrections, 
or higher order l/rric terms, so the total widths do not 



provide a useful bound on Uq^s- 

It is often stated that the difference between the B^ 
and B^ semileptonic rates can be used to constrain the 
impact of WA on the extraction of \Vub\ from B X^^v 

decays. However, Fg^ — Fg^ oc ag, while individually 

and Fg^ , which determine \Vub\^ depend on both 
ao and as. We find no evidence that |ao| <^ \as\, so the 

(B~^^ ( B^^ 

r^L -r^L width difference will not strongly constrain 
the WA contribution to \ Vub\- While the uncertainties in 
our analysis are substantial, it gives strong indication 
that the WA contribution to the B Xu^v rate is less 
than the ~ 3% estimate [|j often used. Our conclusions 
are unchanged if SU{i) breaking or higher order l/mg 
corrections are included, since these will only shift the 
estimate of the WA contribution by ^ 20% of its value. 
The 2a discrepancy in Vut mentioned in the introduction 
cannot be explained away using WA. 

Our results imply that the WA contribution to B de- 
cays, which is a factor {rric/mb)^ ^ 0.03 smaller than 
the corresponding contribution to D decays, is around 
1%. If we use heavy quark symmetry for the bag pa- 
rameters instead of the matrix elements, scaling with 
(r7i^/TO^)(/g//|,) gives 2.5%, still smaller than past es- 
timates. A recent CDF measurement ||3^ of the B me- 
son and Ab baryon lifetimes also indicates that spectator 
effects in the b hadron decays may be smaller than pre- 
viously thought. 
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